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Two recent measurements by the DO Collaboration of spin correlation in tt production using 5.4 fb-l of Teva- 
tron pp collider data are presented. Both rely on the dilepton final state of tt. One measurement relies on full 
reconstruction of the top quark kinematics, and the other makes use of leading-order matrix elements to char- 
acterize the kinematics. The latter measurement is the first ever to have sufficient analyzing power to exclude 
the no-correlation hypothesis. 



1. Introduction 

In the pp — > qq process, the quarks that are produced are unpolarized, but their spins are correlated. This 
is required by angular momentum conservation in the strong interaction. In general, this correlation is unob- 
servable, as the hadronization process involves the emission of gluons that can flip the spins of the quarks. But 
the top quark provides a laboratory for studying the correlation. The short lifetime of top, about 5 x 10^^^ s, 
is shorter than the timescale for strong processes, so top decays before fragmentation and spin flips can oc- 
cur. Thus, the original spin orientation is preserved, and is passed to the decay products. It should then be 
observable through a study of the kinematics of the decay products. 

A measurement of top-quark spin correlation is a test of top-quark properties and also a probe of new physics. 
The very observation of the correlation could in principle be used to set an upper limit on the top lifetime. 
Should top have a non-standard decay (such as t ^ H^b), or non-standard production mechanism (through 
decays of stop pairs, or a Z' resonance), a non-standard correlation would be observed. The correlation is 
ultimately a subtle effect - the theme of this presentation - but there is now enough Tevatron data to explore 
it. 



2. About correlation 

At the Tevatron, the primary production mode of tt is through qq — ti with an s-channel gluon. (This is 
in contrast to the LHC, where the initial state is primarily gg). The q and q must have opposite helicity to 
couple to that gluon, and that forces the t and i to have their spins pointing along the beamline. A correlation 
strength can be deflned based on the number of tt pairs with their spins pointing in the same direction, 

^ ^ Nn + - - ^it ^^ 
^TT + + + ^it ' 

But the spin orientation must be defined with respect to a quantization axis. In the measurements described 
here, the beamline axis, deflned as the direction of the colliding hadrons in the zero-momentum frame of the tt 
system, is used. This choice is intuitive, easy to construct, and optimal for tt produced at threshold. With this 
choice of quantization axis, next-to-leading order QCD calculations predict A = 0.7771q'q42 [I]. 

The spin orientation of the top is then passed to its decay products. The differential angular decay distribu- 
tions of the top daughters is given by 

I dV I, 
r d cos Oi 2 

where cos 9i is the angle between the ith top daughter and the spin of the top. Different decay products have 
different correlation strengths, as indicated by a^. This is illustrated in Figure [l] In the case of a leptonic W 
decay, the lepton has the greatest analyzing power, and for a hadronic W decay it is the down- type quark^ . In 
both cases, a ~ 1. 



In the original observation of parity violation, it was the angular distribution of electrons from nuclear beta decay that were 
observed. What would the history of our understanding of the weak interaction have been had the correlation with the nuclear 
spin direction not been so strong? 
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Figure 1: Illustration of the top decay angles (left) and dependence of the decay rate on the angles (right) [2]. 

Thus, the doubly-differential cross section as a function of the decay angles of decay products from two 
different quarks in tt is given by 



1 



a d COS did COS 02 4 



(1 — Aaia2 cos 9i cos 62)- 



(3) 



To study the spin correlation, one looks for correlation between the directions of decay products from the two 
different top decays. In the case of both tops decaying leptonically, aiai ~ 1, and we write Aaia2 = C ~ A. 



3. Experimental situation^ 

The measurements described below were performed at the Fermilab Tevatron, a pp collider operating at 
y/s ~ 1.96 TeV. Run II of the collider has been in progress since 2001, with almost 12 fb~^ of integrated 
luminosity delivered. The spin-correlation measurements make use of 5.4 fb~^. The data was recorded by the 
DO detector, which consists of silicon and fiber trackers inside a 2 T solenoid, a liquid argon-uranium calorimeter, 
and muon trackers and scintillators inside toroids. 

At the Tevatron, 85% of tt production arises from qq annihilation and the remaining 15% from gg interactions. 
Each top quark decays to Wb nearly 100% of the time, and the final states are characterized by the W decay 
modes. The three final states are all-hadronic, lepton plus jets and dilepton. In the order listed, the final states 
have decreasing rate and increasing number of neutrinos and signal purity. 

The dilepton state in particular is characterized by two high-p^ leptons, missing momentum due to the 
escaping neutrinos, two hadronic jets from b decays, and perhaps additional jets due to initial- and final- 
state radiation. For the purpose of the spin-correlation measurements, the dilepton state has the best analyzing 
power and most accurate measurement of the decay-product (lepton) directions, but the worst statistical power, 
compared to other final states that could be considered. 



4. Event selection 

The two measurements described have the same event selection. A tf-enriched sample is chosen by selecting 
events with two high-p^, isolated, opposite-charge leptons. Only electrons and muons are considered as leptons, 
so the dilepton pairs can be either ee, e/i or fifj,.. There must also be at least two high-p^ jets. To suppress 
backgrounds, a large scalar sum of the lepton and jet px values is required in the e/Lt channel, and significant 
missing energy in the ee and /i/i channels. Backgrounds from Z/"f* (diboson) events are modeled by leading- 
order Monte Carlo samples, and normalized to next-to-next-to-leading (next-to-leading) order cross sections. 



■^This presentation was given at the very end of a long day of top-physics talks, and by unanimous consent of the audience most of 
the material in this section was not presented, as it was deemed redunant. 
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Instrumental backgrounds arise from misidentified tt" and t] decays in electron samples and real muons in jets 
that appear to be isolated in muon samples; both of these are modeled with complementary data samples. The 
selected data sample is about 70% pure in tt events, as shown in Table |lj 

Table I: Estimates of contributions of various physics processes to the selected dilepton sample. 
tt Z/'y* Diboson Instrumental Total Observed 
341 ± 30 93 ± 15 19 ± 3 28 ± 5 481 ± 39 485 



5. Analysis I: Template-based 

DO has performed two different measurements of the spin correlation with this event sample. The first 
is template-based, in which the decay angles are calculated in each event, and the distribution of angles is 
modeled by a sum of templates representing the distributions for correlated and uncorrelated spins [3]. This 
technique has been used before [4], but with a much smaller data sample. 

To observe the correlation, one must measure the angle between the lepton and the beamlinc (which is used as 
the top spin quantization axis) in the zero-momentum frame of the ti system, which requires a full reconstruction 
of the decay. A total of eighteen quantities are needed to specify the final-state configuration, but because of 
the two undetected neutrinos only twelve are measured. Constraining the decay kinematics to the values of the 
top and W masses provides four additional pieces of information, but that still leaves two missing. 

The "neutrino weighting" technique is used to solve the remaining kinematics. Two values of the neutrino 
77 (where 77 = — lntan(0/2) and 6 is the polar angle measured from the beamline) are randomly sampled from 
the neutrino rj distribution as predicted from ti Monte Carlo simulations. These values are then used to solve 
for the implied tt kinematics. This allows a determination of the product of decay angles cos 9i cos 62 and of 
the neutrino momenta. The cos 9i cos 62 value is then weighted by the consistency of the determined neutrino 
momenta with the measured missing transverse energy in the event. The sampling is repeated many times, and 
weighted mean of all solutions obtained is then used as the estimator of cos 9i cos 02 ■ 

ti events are simulated using the MC@NLO generator [B], in which the spin correlation can be turned on or 
off straightforwardly. Then, with the appropriate weighting of simulated samples, cos^i cos02 distributions for 
any value of C can be generated. The left panel of Figure [2] shows the expected distribution of cos 9i cos O2 
at parton level for ti with no spin correlation (C = 0) and standard-model (SM) spin correlation (C = 0.78). 
The distribution is symmetric when there is no correlation, and shifted slightly towards negative values of 
cos 61 cos 6*2 when there is correlation. 




Figure 2: Left: The distribution in cos Si cos 02 for a sample including the NLO QCD spin correlation (C = 0.78) (red 
dashed line) and with no spin correlation (C = 0) (blue solid line) at the parton level, generated using MC@NLO. Right: 
The distribution in cos 6\ cos 62 for the entire dilepton event sample. The summed signal, including NLO QCD spin 
correlation (C = 0.78) (red) and all backgrounds (blue) are compared to data. The open histogram is the prediction 
without spin correlation (C = 0). 

The right panel of Figure [2] shows the distribution observed in the data, along with the appropriately- 
normalized background distribution and the distributions expected for the cases of no correlation and SM 
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correlation. On the face of it, it seems hard to distinguish the two cases. To make a quantitative statement, 
the most hkely value of C is obtained with a binned maximum-likelihood fit. Systematic uncertainties are 
incorporated to the fit as nuisance parameters, and the overall tt cross section is treated as a free parameter to 
avoid biases. 

A Feldman-Cousins-based frequentist approach [5^ is used to set confidence limits on C as a function of the 
measured value. These are shown in FigurejS) The measured value is Cmeas — 0.10±0.45, which can be compared 
with the expected value of C = 0.78. We find —0.66 < C < 0.81 at 95% confidence level. Uncertainties on the 
central value are shown in Table [Til Statistical uncertainties dominate by far; the leading systematic uncertainty 
arises from the limited Monte Carlo statistics in generating the templates. This can obviously be remedied 
when statistical uncertainties are reduced with more data. 




Figure 3: The 68% (inner), 95% (middle), and 99% (outer) C.L. bands of C as a function of Cmcas from likelihood fits 
to MC events for all channels combined. The yellow line indicates the most probable value of C as a function of Cmcas, 
and represents the calibration of the method. The vertical dashed black line depicts the measured value Cmcas = 0.10. 
The horizontal band indicates the NLO QCD prediction of C = 0.777to o42. 



While the measured value of C agrees with the predicted value within two standard deviations, it is also 
consistent with no spin correlation at all. A more powerful technique is required to observe the effect with the 
data sample in hand. 



6. Analysis II: Matrix-element-based 

The second measurement makes use of the leading-order matrix element for tt production and decay, using 
the full event kinematics to determine the fraction of ti events in the sample that the spin correlation that is 
expected in the SM [7]. Matrix elements have never been used in spin-correlation measurements before, and 
their implementation leads to a significant improvement in sensitivity. 

On an event by event basis, we can characterize whether the kinematics are consistent with the SM or with 
no correlation at all. This is done by calculating a probability for consistency of the event with spin correlation 
or non-correlation. The probability is given by 

Psgnix;H) =-^y" /pdf('?i) /pDF('Z2)dgidg2 

. a*, (4) 

Here, Cobs denotes the leading order cross section including selection efficiency, qi and q2 the energy fraction 
of the incoming quarks from the proton and antiproton, respectively, /pdf the parton distribution function, s 
the center-of-mass energy squared of the pp system and d<f>6 the infinitesimal volume element of the 6-body 
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Table II: Summary of uncertainties on Ci 



Source 


+ 


SD 


-SD 


Muon identification 





01 


-0 


01 


Electron identification and smearing 





01 


-0 


01 


PDF 





02 


-0 


01 


Top Mass 





01 


-0 


01 


Triggers 





02 


-0 


02 


Opposite charge requirement 





00 


-0 


00 


Jet energy scale 





01 


-0 


01 


Jet reconstruction and identification 





06 


-0 


06 


Normalization 





02 


-0 


02 


Monte Carlo statistics 





02 


-0 


02 


Instrumental background 





00 


-0 


00 


Background Model for Spin 





03 


-0 


04 


Luminosity 





03 


-0 


03 


Other 





01 


-0 


01 


Template statistics for template fits 





07 


-0 


07 



Total systematic uncertainty 


0.11 -0.11 


Statistical uncertainty 


0.38 -0.40 



phase space. The detector resolution is taken into account through a transfer function W{x,y) that describes 
the probability of a partonic final state y to be measured as x = (pi, . . . ,p„), where pi denotes the measured 
four-momenta of the final state particles. H represents the correlation hypothesis - H — c corresponds to SM 
correlation and H = u corresponds to no correlation. The appropriate matrix element is used for each case. In 
contrast to the template measurement, the full event kinematics plus theoretical models of tt production and 
decay are used, not just the lepton angles, and by adding this information, the sensitivity of the measurement 
is increased. 

For each event, we compute 

Psg„(i? = W)+Psg„(i/ = C) ' 

Events more consistent with having SM spin correlation will tend to have R close to one, while those less 
consistent will have R closer to zero. A value of i? ~ 0.5 implies that it is is difficult to tell which hypothesis is 
more likely. The left panel of Figure [4] shows the expected distribution of R for MC@NLO tt events generated 
with and without spin correlation. In fact, R ~ 0.5 is quite common; the correlation is still a small effect. But 
there is some separation between the two distributions. 

The right panel of Figure [3] shows the distribution in R observed in the data, along with the distribution 
expected for the background events and those for tt with and without spin correlation. By eye, one can see that 
the data are more consistent with the hypothesis of correlation. 

The fraction of events that are consistent with SM spin correlation, /, is obtained from a binned likelihood fit 
that is very similar to that of Analysis I. This fraction is of course expected to be 100%. The confidence bands 
are shown in Figure [H] We find / = 0.74^0 45, which is consistent with the SM. A value of / = is excluded 
at the 97.7% confidence level. (From ensemble testing, the measurement was expected to exclude / = at 
the 99.6% confidence level; if the correlation exists as expected in the SM, this measurement was " unlu cky" 
in finding a value less than 100%.) The uncertainties on the measured value of / are listed in Table III once 
again, statistical uncertainties greatly dominate. 



7. Summary 



Quark spin correlation is a phenomenon that can only be seen in tt production, thanks to the short top 
lifetime. However, it is a subtle effect that requires large data samples and sophisticated analysis techniques to 
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Figure 4: Left: Comparison of the discriminant R between SM spin correlation H = c and no spin correlation H — u at 
parton level. The first and last bin include also the contributions from R < 0.29 and R > 0.63. Right: The predicted 
discriminant distribution R for the combined dilepton event sample for the fitted au and /meas compared to the data. 
The prediction with spin correlation (/ — 1) and without spin correlation (/ = 0) is shown including background. 




meas 



Figure 5: The 68.0% (inner), 95.0% (central), and 99.7% (outer) C.L. bands of / as a function of /meas from likelihood 
fits to MC events. The thin yellow line indicates the most probable value of / as a function of /meas, and therefore 
represents the calibration of the method. The vertical dashed black line indicates the measured value /meas = 0.74. 

observe. Indeed, the matrix-element technology is perhaps the most powerful, and most complex, analysis tool 
that is currently available for Tevatron data analyses, and it was required here to have the hope of observing 
the effects of interest. Two analyses of dilepton events at DO have been performed. One was a template-based 
analysis using full reconstruction of top decays, giving a result within two standard deviations of the NLO QCD 
prediction, but also compatible with the no-correlation hypothesis. The other one was a matrix-element-based 
analysis that gives a result consistent with the SM hypothesis, and powerful enough to exclude the no-correlation 
hypothesis for the first time ever. Both analyses are statistics limited, with only about half of the final DO Run II 
data sample analyzed so far. Thus, there is great potential for improving the precision of the measurements in 
the near future. 
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Table III: Summary of uncertainties on /, 


me as • 


Source 


+SD 


-SD 


Muon identification 


0.01 


-0.01 


Electron identification and smearing 


0.02 


-0.02 


PDF 


0.06 


-0.05 


mt 


0.04 


-0.06 


Triggers 


0.02 


-0.02 


Opposite charge selection 


0.01 


-0.01 


Jet energy scale 


0.01 


-0.04 


Jet reconstruction and identification 


0.02 


-0.06 


Background normalization 


0.07 


-0.08 


MC statistics 


0.03 


-0.03 


Instrumental background 


0.01 


-0.01 


Integrated luminosity 


0.04 


-0.04 


Other 


0.02 


-0.02 


MC statistics for template fits 


0.10 


-0.10 


Total systematic uncertainty 


0.15 


-0.18 


Statistical uncertainty 


0.33 


-0.35 
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